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Abstract 

The twist three contributions to the Q 2 -evolution of the spin-dependent structure function 
92(xBj,Q 2 ) are considered in the non-local operator product approach. Defining appropri- 
ate twist three distribution function we derive their evolution equation for the nonsinglet 
case in leading order approximation. In the limit XBj —> 1 as well as in the large N c 
limit we confirm the result that the evolution of the nonsinglet part of 52 is governed by a 
Gribov-Lipatov-Altarelli-Parisi equation. 



1 An earlier version of this paper has been presented at the "3rd Meeting on the Prospects of Nucleon-Nucleon 
Spin Physics at HERA", JINR Dubna, 28-29.6.1996. 



1 Introduction 



Recently, in deep inelastic scattering the first moments of the polarized structure function 
g 2 (xBj,Q 2 ) are measured |]. In leading order of the momentum transfer Q 2 this structure 
function is determined by twist two as well as twist three contributions. In comparison with 
the twist two case the leading order analysis for the twist three part is more subtle due to 
the appearance of a set of operators mixing under renormalization and constrained by relations 
between themselves. 

Up to now there exist already several papers which determine the local anomalous dimensions 
^ |J and the evolution kernels for the distribution functions as well as nonlocal operators 
0, |5], |6|. In the nonsinglet sector the one- loop result || for the evolution kernel of light-ray 
operators was confirmed in 0. It has been checked that the local anomalous dimensions coincide 
with the results given in || . 

The renormalization properties of twist three operators indicate that the evolution equa- 
tion for the distribution function is more complicated than the Gribov-Lipatov-Altarelli-Parisi 
(GLAP) equation govering the evolution of twist two parton distribution functions. However, for 
the nonsinglet case the solution of this equation is known in both limits N c —>■ oo and Xbj —> 1 
0- 

Here we introduce a new twist three distribution function defined in terms of three-particle 
operators and we give first results about the evolution in the flavor nonsinglet sector. This 
distribution function depends on an effective momentum fraction and on the position of the 
gluon field. The evolution is governed by an extended GLAP equation. The moments are 
derived from twist three operators with definite spin and so that they do not mix with each 
other. The remaining mixing problem can be treated numerically for the first few moments. In 
the limit N c — > oo this equation reduces to an evolution equation of the GLAP type. 

2 Nonlocal operator product expansion 

The antisymmetric part (with respect to Lorentz indices) of the hadronic tensor 
W* = e^xsq 

being relevant for polarized deep inelastic scattering is given by the structure functions 
gi(xBj, Q 2 ) and giixBji Q 2 ) depending on the Bjorken variable xsj and the momentum transfer 
q, q 2 = —Q 2 . The polarization vector of the proton is defined as S s = u(P, S)'j S 'y 5 u(P, S). 
Our conventions are 7 5 = Z7°7 1 7 2 7 3 and 60123 = 1- The hadronic tensor is determined by the 
absorptive part W^ v of the virtual forward Compton amplitude T^ v , 

T^ = - [ d 4 xe ixq (P, S\Tf(x/2)f(-x/2)\P, S). (2) 

7T J 

In the Bjorken region the leading terms in Q 2 of the Compton amplitude correspond to the 
leading light-cone singularities in the coordinate space so that the light-cone expansion for op- 
erator products (OPE) can be applied. This provides the factorization of the structure function 
in a perturbatively determined coefficient function and a nonperturbative parton distribution 
function. 



^pgi(x Bj , Q 2 ) + qPS f qP y P 92(x Bj , Q 2 ) 
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The leading order analysis starts with a perturbative investigation of the time ordered product 
of two electromagnetic currents. (In the following we neglect flavor and color indices): 

{Tf(x)f(y)r = e^i&" p D c (x, y) $(x)w 5 U(x, y)^(y) + (x, fi <- yu) + ■ ■ ■ , (3) 

where the path ordered phase factor U(x,y) = Pexp j— ig J* d«A4 At («;)} ensures gauge invari- 
ance. For simplicity we set y = 0. 



Technically, we apply the light-cone expansion [11, O] proved by Anikin and Zavialov in 



renormalized quantum field theory. Heuristically, this nonlocal OPE is obtained by approximat- 
ing the vector x by the light-like vector x, defined as x = x + a(x,rj)rj, where rj denotes a fixed 
auxiliary vector (e.g. normalized by rf = 1) and a the corresponding coefficient. In leading order 
we substitute x —* x whereas the ^-singularities remain unchanged, i.e., x 2 — * x 2 . In general, 
the coefficient function will depend on two auxiliary variables whose range (according to the 
a-representation of the contributing Feynman diagrams) is restricted by < Kj < 1. Here we in- 
troduce these variables quite trivially through integration over two ^-functions. In this intuitive 
way we get the following light-cone expansion: 

{Tf(x)f(0)r = j^dm ^d^^-m^d^D'^iO^K^ + iK,^^)},^) 

O p (k 1 ,K 2 ) = i)(K 1 x)j^ 5 U(KiX,K 2 x)'4)(K2x). (5) 

The light-ray operator O p contains both twist-2 and twist-3 contributions. In the following 
we use the light-cone gauge xA = so that U(k\X, k 2 x) = 1. The twist-2 part of O p (ki, K 2 ) can 
be obtained by contraction with the light-cone vector x: 

O tw2 (/ti, k 2 ) = x p O p (k u k 2 ) = ijjfax) $j 5 ip(K 2 x). (6) 

It turns out (for more details see 0) that for an application to forward scattering the twist-3 
part can be written as 



O t p w3 (Ki,K 2 ) = —i(K 2 — Ki) / duuO p (Ki, K\U + K 2 U) 

J 



O p {kx,k 2 ) = I duipiKx^i [j p x a - ^5f p(7 ]7 5 L' fT (u, Ktx, k 2 x)i/;(k 2 x) , (7) 







where we introduced the notation D p (u, KiX, k 2 x) = d p K2i . + igA p ([niu + H2u]x) and u = 1 — u. 

Using the results of the OPE (§-|7|) it is straightforward to obtain the structure function in 
terms of parton distribution functions. As is well known the longitudinal structure function g\ 
reads 

9i(x Bj ,Q 2 ) = i J2 e 2 q Aq q (x Bj ,Q 2 ). (8) 

q=u,d,.. 

The quark distribution functions (containing quark and antiquark contributions) are defined as 
matrix elements of leading twist-2 operators, 

Aq q (x, Q 2 ) = J (P, S |o tw2 (0, + -k) j P, S) e™« sp \ (9) 

where the renormalization point of the operator is set equal to the momentum transfer Q 2 . The 
leading order analysis ||,|| suggests that the structure function g 2 can be written as 

92(x, Q 2 ) = Q 2 ) + T - Uv, Q 2 ) (io) 

Jx y 
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where g 2 (x,Q 2 ) = gi(x,Q 2 ) + g 2 (x,Q 2 ) is decomposed into the twist-2 part given by g\ and a 
remaining twist-3 part g 2 . The twist-3 contribution 

U*B V Q 2 ) = \ E e 2 q Aq q (x Bj ,Q 2 ), (11) 

q=u,d,.. 

can be expressed by the distribution function 

P, S) e lXKi * p \ (12) 
where, as will be seen below, unlike q q (x,Q 2 ) this function has no simple parton interpretation. 



Aq q (x,Q 2 



S P (RS 



3 Evolution kernels for twist-3 light-ray operators 



The twist-3 operator O p is not closed under renormalization. Indeed it will mix with three- 
particle operators which contain also the gluon field strength. Using the equation of motion, 
(i p — m)ip = the operator O p can be decomposed with respect to gauge invariant twist-3 
operators: 



O p (k!, K 2 ) 



[K 2 ~ Kl 



du{ — 2M p (ki, k\u + k 2 u) + u + S p (k,i, k\u + k 2 u, k 2 ) 

+ u ~S p (ki, K\U + K 2 U, K 2 )} H . (13) 



Here ^Sp are the nonlocal generalizations of the so-called Shuryak-Vainshtein operators O 



iz S p (K 1 ,T,K 2 ) = igip{Kix) $ iF ap (rx) ±-f 5 F ap (rx) x a i){n 2 x) 



where F, 



i 



a/3 2 ^a/3fiu 



F^ u is the dual field strength tensor. Furthermore, 



M p (ki, k 2 ) = m ^(Kix)a ap x a ^(xD)(K 2 x)tp(K 2 x), a afi = ^[y a , jp\ 



(14) 



(15) 



denotes a mass dependent operator. Besides the equation of motion operators we neglected 
also trace terms (proportional to x p ) and operators which vanish in the forward case (from 
general principles |TJJ] it is known that the equation of motion operators do not contribute to 
the evolution of physical matrix elements of gauge invariant operators). 

For simplicity we will neglect the mass operator. Then the operators S p which are related to 
each other by charge conjugation are closed and, furthermore do not mix under renormalization. 
We calculated the evolution kernels for these operators in light-cone gauge using the Leibbrandt- 
Mandelstam prescription | 15| . Since both kernels are related to each other we give only the result 

d , a s (fi 2 ) 



for S p : 



dfi 2 



S P (k!, K 2 ) 



4tt 



dy J V dz{(2C t 



y5{z) S p {ki - n 2 y, -K%y) 



+C/ 



K{y,z) 



1 + S(z) V - + S(y) - 
y z 



-2z S p (-K 1 y, k 2 - K\z) + K(y, z) S p (K 1 y + K 2 y, k 2 z + K\z) 
(2z + L(y, z)Ys p (n 1 y, k 2 - k^z) + L(y,zyS p (ni - K 2 z,K 2 y) }, (16) 
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L(y,z) 



5(1 - y 



z) — + o(z) - 

y 1 v 
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S(y)S(z). 



To condense the notation we used ± S p (ki, k 2 ) = ± S p (ni, 0, k 2 ) and the standard plus-prescription 
fulfilling Jdy[...} + = and J dydz[...] + = 0, respectively; C F = (N 2 - 1)/(2N C ) = 4/3 and 
Ca = N c = 3 are the usual Casimir operators of SU C (3). 
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4 Definition of twist-3 parton distribution functions 



Since the twist-3 operators are in fact three-particle operators it is also necessary to extend the 
definition of the distribution function ([12]) . A quite natural way is to take the Fourier transform 
of the operators ± S p (ki,k 2 ) sandwiched between polarized proton states. The resulting C-even 
distribution function 

AQ NS ( Xl ,x 2 ) = J ^^^^ 2 ( P ,S\ + S p (k 1 ,k 2 ) + -S p (k 2 ,k 1 ) + {k 1 ^k 2 }\P,S) (17) 



v IKIXX {xP)+iK2X2 ixP) 



depends on the momentum fractions x%, x 2 and posseses the support property \xi\ < 1, |xi =bx2 1 < 
1. Unfortunately, the resulting evolution equation will be very complicated 0, ||. To get a 
simpler one and to be able to diagonalize it with respect to one variable we take into account 
that operators with different spin do not mix with each other. 

C-even operators with definite spin can be obtained easily from Shuryak-Vainshtein operators 

by 

Y"{k,u) = + S p (0,ku,k) + ~S p (k,ku,0) (18) 

k=l \ / 



The operators Y p (u) possess the definite spin n. In Eq. (|IBD the variable u gives (for k — 1) the 
position of the gluon field on the light cone. For < u < 1 the gluon field lies between the two 
quark fields. Note that the gluon field can also be outside of this range so that u is not restricted 
to the interval [0,1]. 

Consequently we define a new distribution function as Fourier transform with respect to the 
variable k only, so that 

A<f s (y,u) = J^--^-(p,S\Y p (K,u)+Y p (-K,u)\P,S)e i ^ (20) 

depends on the Fourier conjugate variable y and the gluon position u. To clarify the meaning 
of the variable y we express this distribution function ( p0|) in terms of the distribution function 

AQ NS (x u x 2 ), 1 

Aq NS (y,u) = J dx 1 J dx 2 9(1 - \x 1 ± x 2 \)S(y + ux 1 - mx 2 )AQ ns (xi, x 2 ) (21) 

which gives the support restriction \y\ < Max(l, \2u — 1|), i.e., for < u < 1 the variable y can 
be interpreted as an effective momentum fraction. 

The evolution equation for the distribution function (|2H|) can be derived from the renormaliza- 
tion group equation (|T6|). Using the transformation (|T8|) and the definition (|20"D a straightforward 
calculation provides the following result: 

g2 ^ Ag ~ NS ( y ' M ' g2 ) = ^^S d il dvp{z ^ v) ^^ Q2 )- (22) 

Here the integration region is determined by both the support of Ag NS and the support of the 
evolution kernel 
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Pqq(z,U,v) 

2Cf — Ca 



(23) 



[@i(z, u, v)K(z, u, v)] + + 6 2 (2, u, v)L(z, it, v) — ©2(-2, u, v)M(z, u, v) 



+ 



Ca 
2 



[Q 3 (z,u,v)N(z,u,v)] + + 



it — > It 



+ 03(2, ix, d)M(z, rx, d) — -<5(ix — d — -2) 



where [A (2, ix, u)] + = it, 1;)— 5(1— z)5(u— v) Jq dz' J dv' A(z',u,v') and the auxiliary functions 
are defined by: 

Qi(z, ix, v) = 9(z)9(u — zv)9{u — zv), Q2(z,u,v) = 9(—uvz)9({l — vz}u)9({z — u}u), 
Q 3 (z,u,v) = 9(uz)9(uvz)9({vz — u}u), 



K(z, ix, v) 
M(z,u,v) 



u cv x / It — > It 

Z + < — rOflX — ZV) + 

DID — IX) \ V — > V , 



vu 2 



, L(z,u,v) = -sign(u)— T 5(u - z) 



2z(l - zv) 



U" 



at/ x sign(ix)D I v , . u . , . 
, N(z, ix, v) = \ < -6(1 -z) + —5(u- zv) 

u(v — ix) I IX V 



(24) 



5 Solution of the evolution equation 

A diagonalization of the evolution equation fl22|) can be achieved partly by the Mellin transform 

Agf , u \ = f 1 dx x n-2 A ~NS , ^ _ (25) 
JO 

These moments are given by matrix elements of the operators Y^(u) with spin n and are poly- 
nomials in the variable u of order (n — 2) [see Eq. (pLQf) 1 . Their evolution equation reads 

Q 2 ^Agf(ix,Q 2 ) = ^^fdvP n (u,v)Aq™(v,Q 2 ), (26) 

where the n-dependent kernel P n [u,v) = Jq dz z n ~ 2 P n (z,u,v) can be calculated from Eq. (|23|) . 
For each given n, Aq^ s (u) can be expanded with respect to the (n — 2) eigenfunctions e n> j(it) of 
the kernel P n (u,v). These polynomials can be constructed by diagonalization of 

("_"i)/ d ^;(^)^ H = e("I 1 2 )^« b ^- 1 « < - 1 . (27) 

where 7™- are the known anomalous dimensions of the local operators [H HI- A general 
analytical solution of this problem is not known. But, for the first few moments the numerical 
solution can be obtained easily. 

Knowing the [n — 2) eigenfunctions e n> i(u) and eigenvalues X nii of the kernel P n (u,v) it is 
possible to determine the evolution of the moments 

Agf(ix,Q 2 ) = |> (Q 2 ) e„>) exp ( y^A n l (28) 

i=l [ J Qo ^ Z7t ) 
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The nonperturbative coefficients q (Qq) at a reference momentum squared Qq have to be deter- 
mined from the initial value Ag,^ s (u, Qq). Since, assuming the validity of Eq. (PH|), the twist-3 
nonsinglet part of the structure function g 2 satisfies 

x~g^ s (x) = ^Aq NS (x), Aq NS {x,Q 2 ) = f 1 duuAq NS {x,u,Q 2 ) (29) 
dx Jo 

it is obvious that the initial values can not be obtained alone from transversal polarized deep 
inelastic scattering experiments. 

Finally, we show that the above mentioned complications do not appear in the limit x — > 1 as 
well as in limit of large Nq. We deal directly with Aq NS (x). From the definition of the evolution 
kernel ( p3|) it turns out that J X duuP(z, u, v) factorizes in boths limits according to vP(z). Thus, 
from the evolution equation ( [22] ) follows for x — > 1, 

(30) 

Ca) 5(l-^)|Ag NS (^,Q^+ (l-x), 
and for large N c , 




These equations are of the GLAP type and coincide with the result of flOfl . The 1/N C suppressed 
terms contain contributions from the region v < and v > 1 and they provide a system of 
coupled evolution equations. 
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